Particle number fractionization of an atomic Fermi-Dirac gas in an optical lattice 
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We show that a dilute 2-species gas of Fermi-Dirac alkah- 
metal atoms in a periodic optical lattice may exhibit frac- 
tionization of particle number when the two components are 
coupled via a coherent electromagnetic field with a topologi- 
cally nontrivial phase profile. This results in fractional eigen- 
values of the spin operator with vanishing fluctuations. The 
fractional part can be accurately controlled by modifying the 
effective detuning of the electromagnetic field. 

PACS: Q3.75.Fi,Q5.30.Jp,11.27.-Hd 



Particle number fractionization is a remarkable phe- 
nomenon in both relativistic quantum field theory and 
condensed matter systems Jackiw and Rebbi [||^] 

showed that for a fermionic field coupled to a bosonic 
field with a topologically nontrivial soliton profile, the 
fermion number can be fractional. The noninteger par- 
ticle number eigenvalues may be understood in terms of 
the deformations of the Dirac sea (or the hole sea) due 
to its interaction with the topologically nontrivial envi- 
ronment. In this paper we propose a manifestation of 
this phenomenon in the atomic regime, using an optically 
trapped Fermi-Dirac (FD) atomic gas ||]. Fractional 
fermion number has been demonstrated previously in the 
condensed matter regime in ID conjugated polymers by 
Su, Schrieffer and Heeger |^,0,D . Fractionally charged ex- 
citations are also fundamental to the fractional quantum 
Hall effect ||^ , but the fractionization mechanism is very 
different from that in the polymers and in the atomic 
gas in this paper. Our dilute atomic gas has a possi- 
ble advantage, compared to condensed matter systems, 
in the sense that the interatomic interactions are weak 
and well-understood, and there now exists a wide range 
of atomic physics technology to detect, manipulate, and 
control atoms by means of electromagnetic (em) fields. 

We study a two-species atomic FD gas in a ID optical 
lattice, coupled to an em field with topological properties 
similar to a soliton, or a phase kink. We show that this 
topologically nontrivial coherent background field results 
in fractionization of the particle number operator eigen- 
values for the fermionic atoms. Also, the spin operator 
has fractional eigenvalues with vanishing fluctuations. 

The background field is generated by means of a co- 
herent em field inducing transitions between the two 
fermionic components occupying different internal lev- 
els, e.g., by using experimentally realized technology of 
rapidly rotating laser beams 1^,0. In the low energy 
limit we demonstrate the one-to-one correspondence of 



the Hamiltonian for the FD atoms to the Jackiw-Rebbi 
relativistic Dirac Hamiltonian describing fractionization 
in quantum field theory. This is related to fractioniza- 
tion in the polyacetylene polymer systems, where the lin- 
earized lattice vibrations are coupled to the electron dy- 
namics which becomes analogous to that for a relativistic 
Dirac equation exhibiting fractional particle number. We 
also show how our proposed system could be generalized 
to higher spatial dimensions to represent the fractioniza- 
tion in relativistic 2-1-1 and 3-1-1 D quantum field theories. 

Tremendous progress in experiments on cold trapped 
alkali-metal atomic gases has allowed numerous studies 
with Bose-Einstein condensates and the cooling of a FD 
gas to the quantum degenerate regime jl^. Recently, 
condensates have been loaded to a periodic optical lattice 
, a system which is expected to exhibit interesting dy- 
namical phenomena [T^ ]. Similar experimental progress 
is anticipated for FD gases. 

The vacuum state in an atomic condensate can support 
topological excitations including defects, such as solitons 
[1 and vortices [Q, as well as SU(2) Skyrmion textures 
In this paper we show that in atomic gases also 
the vacuum state itself may display nontrivial topological 
quantum numbers with a strong analogy to the vacua 
encountered in relativistic quantum field theories. 

The phenomenon of fractional fermion number is best 
illustrated by the following 1-1-1 dimensional Dirac hamil- 
tonian for a 2-component spinor ^/{x) coupled to a 
bosonic condensate (p{x), which can be taken to be a 
static classical background field: 



H 



dx [c?i*V^-^ -f hgip-^'^a^-^ + mc^*V^*] . (1) 
idx 



Here cr* denote the Pauli spin matrices, g the coupling 
coefficient, m the fermionic mass, and c the velocity of 
light. We show below that the atomic FD hamiltonian 
can be written in this form. 

We assume the bosonic field has a doubly-degenerate 
ground state with constant field ^p{x) = ±7. The vacuum 
then exhibits a spontaneously broken refiection symme- 
try if <-> —if. For m = 0, the Dirac Hamiltonian has 
a charge conjugation symmetry, so that for every eigen- 
value e there exists an eigenvalue — e, and the correspond- 
ing eigenfunctions are paired according to = cr^^*. 
The fermion particle number operator is 



1 



dx[^\x),^{x)] 



(2) 



In the free vacuum, the fermion particle number vanishes. 
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The soliton background (p{x) interpolates between the 
two vacua: (p{co) = — </?(— oo) = 7. For m = 0, charge 
conjugation symmetry is preserved, so that positive and 
negative energy states are paired, but in addition to con- 
tinuum modes there is now also a zero-energy bound state 
localized at the soliton jump. This state is charge self- 
conjugate and results in a doubly-degenerate soliton sec- 
tor vacuum. The number operator in the presence of the 
soliton reads |^: 

N = a^a- 1/2 + J dk {blbk - 4cfc) , (3) 

where bk and Ck denote annihilation operators for contin- 
uum fermion and antifermion modes (respectively) , while 
the operators a and couple the two degenerate zero- 
energy ground states. The ground-state soliton states 
possess fractional fermion numbers ±1/2. The Hamil- 
tonian is diagonal in the number representation and all 
the fermion eigenstates display half integral eigenvalues 
with vanishing fluctuations. The fractional part of the 
fermion number has a topological character: it is insensi- 
tive to local deformations of the bosonic field, depending 
only on its asymptotic behavior. For m 7^ 0, charge con- 
jugation symmetry of the Dirac Hamiltonian (|^) is bro- 
ken and the positive and negative energy states are no 
longer coupled in a simple way. The particle number of 
the soliton vacuum is (0|Af|0) — — l/7rarctan (Tig^/mc^) 
and may exhibit arbitrary fractional eigenvalues [|l5| . 

In our scheme to realize particle number fractionization 
in atomic gases we consider neutral FD atoms loaded in a 
periodic optical lattice. The confining optical potential is 
induced by means of the ac Stark effect of the off-resonant 
laser beams [|2| . We assume a FD gas with two internal 
levels t and i coupled via an em-induced transition. The 
coupling could be a far-detuned optical Raman transition 
via an intermediate atomic level, a microwave, or a rf 
transition. Furthermore, we assume that the two species 
experience optical potentials which are shifted relative to 
each other by A/4, where A is the wavelength of light of 
the confining optical lattice. This is realized, e.g., when 
the laser beam is blue-detuned from the internal tran- 
sition of the atoms in level "f, and rcd-dctuned by the 
same amount from the internal transition of the atoms 
in level |. A simple example ID lattice potential in that 
case is V^{x) = Vb sin^(fca::), and Vi{x) = — Vbsin^(fcx). 
The neighboring lattice sites represent atoms in different 
internal levels and are separated by a distance A/4. The 
Hamiltonian for this two-species FD gas is 

i k odd 

- E (W+iA+H.c.) . (4) 

/ even 

Here ai and f3j denote annihilation operators for atoms 
in levels t and |, at lattice sites i and j, respectively. The 
em-induced coupling between the two internal states is 



described by K = / d^r^'^{r — rk)i^{r)tpi{r—rk±i), and 5 
stands for the effective detuning between the levels. The 
em-coupled terms are the analogues of the hopping terms 
in the corresponding polymer hamiltonian Q . The mode 
functions of the individual lattice sites (Wannicr func- 
tions) arc denoted by ipj{r — r^). We assume that the em 
coupling between the internal levels with frequency f2(r) 
is the only transition mechanism for the atoms between 
neighboring lattice sites and therefore we ignored the di- 
rect tunneling. For simplicity, we also ignore the s-wave 
scattering between the two FD species. 

To produce fermion number fractionization, we pro- 
pose to take the couphng frequency r2(r) to be a phase- 
coherent superposition of a standing em field along the x 
axis and a field with the spatial profile ip{x): 

n{r) = iV{r) [sin(2fca;) + ip{x)] , (5) 

with k = 27r/A We show that such a coupling fre- 
quency converts the atomic lattice Hamiltonian into 
the Dirac Hamiltonian (|l|) in the continuum limit. This 
is similar to the ID polymer case but the physics 

is very different: in our atomic system the kink (p{x) ap- 
pears in the em-induced coupling between the internal 
atomic states, and not as a physical domain wall kink. 




FIG. 1. The energy diagram of the two-species Fermi-Dirac 
gas in an optical lattice. The atoms occupy two different in- 
ternal levels and experience different periodic optical poten- 
tials shifted by A/4 (on the left). The coupling between the 
neighboring lattice sites, representing different spin compo- 
nents, is induced by two-photon Raman transitions, or by a 
superposition of the Raman transition and a one-photon mi- 
crowave or rf transition (dotted line). The transition region 
is denoted by the dark shaded area representing the overlap 
region of the atomic wave functions. The energy difference 
between the two components in the overlap region is denoted 
by ui. The two-photon transition is far-detuned by A from an 
intermediate atomic level (on the right). 

For simplicity, we assume that ri{x) = J d'^rip^(r — 
rfc)V(r)?/;| (r — rk±i) does not change its sign over the 
length of the lattice. The purpose of the standing wave 
in Eq. (^ is to introduce an alternating sign between the 
neighboring lattice sites. Such a coupling may be pre- 
pared by a two-photon optical Raman transition (Fig. |^). 
The strength of an off-resonant two-photon Rabi fre- 
quency in the limit of large detuning. A, from the in- 
termediate state is f2 oc 7?.i7?,2/A, where TZi denote the 
Rabi frequencies in the individual transitions ||l7| . For 
two standing-wave one-photon couplings displaced from 
one another by A/4, the two-photon Rabi frequency is 
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n cx sin(fca;) cos(fca;) cx sin(2fca;). 

We choose the field profile ip{x) to exhibit a phase jump 
of TT at X 0. The phase jump represents a topological 
phase singularity, or a phase kink. This type of a coupling 
might be produced either by making use of the (rf or mi- 
crowave) transition between the spin states, or by means 
of an optical Raman transition other than the one used to 
produce the sin(2A:a;) standing wave. A phase profile with 
topological properties similar to ip{x) could be prepared, 
e.g., by means of a standing microwave cx sm{qx) with 
q<^k. It should also be possible to shape the wave fronts 
of the coupling lasers to produce a two-photon transition 
with a desired phase jump. To avoid rapid phase vari- 
ation at the length scale A, one could use either laser 
beams co-propagating along the x axis, or beams with 
the wave vectors nearly perpendicular to x. Alterna- 
tively, the coupling could possibly also be obtained by 
a dc magnetic field B{x) cx x as explained in Ref. [ p^ . 
An em field, topologically similar to ^p{x), was also used 
to create solitons in atomic condensates Q. Unlike in 
Ref. 1^ we could assume that the coupling field itself is 
formed by two rapidly rotating em fields resulting in a 
desired time-averaged phase profile. 

We assume that sin(2A:x) and (p{x) in Eq. (||) are ap- 
proximately constant over the spatial overlap area of 
neighboring lattice site atom wavefunctions. Then the 
sine function is approximately ±1 at each overlap area: 

K^iri[{-ir + ^{x)\. (6) 

For notational simplicity, we take /t* — — k, and "q real. 

In this paper we study the Hamiltonian (|^) [with ri(r) 
defined in Eqs. (||) and (^)] only in the continuum limit, 
where it can be transformed to the relativistic Dirac 
Hamiltonian ([^) exhibiting fractional charge It may 
be shown that the fractional eigenvalues also emerge in 
a more general case, but the continuum field theory is 
amenable to a simpler description. The continuum limit 
corresponds to the linearization of the fermionic band 
structure and becomes accurate in the dilute gas limit, 
where the atomic correlation length is much larger than 
the lattice spacing. In the continuum limit we write the 
fermionic annihilation operators as continuous functions 
of the lattice spacing d= A/4: 

a, = V2du{jd), = V2dv{jd) . (7) 

Then the continuum limit proceeds exactly as in the poly- 
mer case [Q. To leading order in small d we obtain 

H/h = 2id^r]^[u''{nd)v'{nd) +v^{nd)u'{nd)] 

n 

Sd 

+ — ^y^Ju^ {nd)u{nd) — v\nd)v{nd)] 

n 

+2irid ^ (p{nd) [u^ {nd)v{nd) -|- {nd)u{nd)] . (8) 

n 

Here u'{nd) denotes a discrete spatial derivative of u. In 
the continuum limit we replace nd ^ x and dJ2n ~* 



J dx. By introducing the spinor ^{x) = [u{x) v{x)]'^ and 
the transformation: 5* — > exp {ina^ / 4)^/ , we may express 
Eq. as the relativistic Dirac Hamiltonian p), when 
we identify c = A77/2, g — 2ri, and m — 2h5 / {)? if') . In 
this case the spinor components refer to the two internal 
atomic levels. Note that for nonzero detuning 5 the sys- 
tem is not charge conjugation symmetric, and the eigen- 
values can have any fractional value. The ratio between 
the coupling strength 777 and 5 therefore determines the 
fractional part of the particle number. In experiments 
this could be engineered accurately, allowing a controlled 
way of preparing the fractional part of the eigenvalues. 

The crucial part of our proposal for fractional particle 
number is the em field (^(x) in Eq. (Q). This is very dif- 
ferent from the fermion particle number fractionization 
in polymers, as our fermionic fields are not coupled to a 
bosonic matter field with a domain-wall soliton. Instead, 
the coherent em field with a topologically appropriate 
phase profile is coupled to the FD atoms via internal 
transitions. This results in the quantization of the FD 
atomic gas with nontrivial topological quantum numbers 
corresponding to the soliton sector of the relativistic 1-1-1 
quantum field theory models of fractionization. On the 
other hand, a spatially constant phase profile ^{x) rep- 
resents the FD vacuum sector exhibiting integer particle 
numbers and no bound state. 

The normalizable bound state, which plays an impor- 
tant role in the fractionization, belongs to only one of the 
two fermionic components, independently of the shape or 
the position of the phase kink: for the solitonlike phase 
kink, f{h) = —ip{—h) — 7, for all h ^ d, only level 
"f is occupied, and for the antisolitonlike phase kink, 
(p{h) = —ip{—h) — —7, only level | is occupied 
When TO = 0, the bound state has zero energy. Its spa- 
tial profile ^ exp (— 7|a;|/d) (for a sharp kink at x = 0) 
depends on the relative strength of the superposed em 
fields, according to Eq. (^), determining 7. Unlike in the 
polymer case, where the size of the bound state is fixed, 
in the atomic case this could be varied experimentally. 

Because the local density of states is conserved, a zero- 
energy mode creates a fractional deficit of states in both 
the valence and the conduction bands. In the presence 
of charge conjugation symmetry, the density of states is 
a symmetric function of the energy, and both bands have 
a decifit of one-half a state. By assigning the atoms in 
the conduction and valence bands as "particles" and "an- 
tiparticles" , respectively, we can interpret the fractional 
particle number operator as the occupation number dif- 
ference between the bands. However, because the zero- 
mode always occupies only one spin component at a time, 
also the spin is fractionized. As an example, the two 
species may correspond to the eigenstates of a single- 
particle spin operator along the z axis, az, with eigen- 
values ±1. Then the eigenvalues of the many-particle 
operator Sz = J2n localized around the soliton, are 
fractional with vanishing fluctuations. 

The total number of atoms, of course, must remain 
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an integer and any realistic optical lattice has a finite 
size. For every fractional particle number located at the 
phase kink (forming a soliton), some fractional charge is 
distributed at the boundary of the atomic cloud, or is 
associated with an antikink. Although the fluctuations 
of the fractional eigenvalues do not then in a finite lat- 
tice exactly vanish, it can be shown that the fluctuations 
decay exponentially as a function of the size of the sys- 
tem and can be considered negligible, if the size is much 
larger than the atomic correlation length There- 
fore, every localized measurement of the particle number 
around the phase kink can yield a fractional result. 

In experiments on fractional fermion number we may 
detect the bound state or measure a fractional expec- 
tation value and determine its fluctuations to ascertain 
that they are compatible with fractionization. The FD 
gas exhibits a gap ^ 2r/ in the excitation spectrum. The 
phase kink creates a bound mode at the center of the 
gap, hence excitations at half the gap energy |^0|. These 
midgap transitions Q could be probed in resonance spec- 
troscopy. The bound state also alters the dynamical 
structure factor, which may be observed via light scatter- 
ing JTzt . The optical signal may be magnified by simul- 
taneously preparing many phase kinks. The fractional 
particle number could be detected by measuring the oc- 
cupation numbers of the individual lattice sites. For in- 
stance, a magnetic field gradient may be introduced that 
causes a detectable change in the spin flip frequency from 
site to site Finally, an off- resonance optical probe 

couples to atom density Fluctuations in the scat- 

tered light should therefore convey information about not 
only the particle number, but also about its fluctuations. 

We can also generalize the proposed scheme to particle 
number fractionization in higher dimensional models in 
2D or 3D optical lattices with atoms coupled to em fields 
exhibitingphasc profiles similar to topological defects or 
textures |14| . In relativistic 2-f 1 D quantum field theory, 
a fermionic field coupled to a bosonic field exhibiting a 
vortex profile results in fermion particle number fraction- 
ization Q. In our approach this would correspond to a 
2D optical lattice of FD atoms coupled by an em field 
with a nonvanishing phase winding around any closed 
loop circulating the axis of a vanishing field amplitude. 
This is also similar to the Raman field used in the JILA 
experiments to couple two atomic condensates in order 
to create vortices ||T^. We may find an analogy to rel- 
ativistic 3-1-1 D quantum field theory in the presence of 
a nonvanishing SU(2) topological charge of the bosonic 
field by means of the em field configurations proposed in 
Ref. [h4 to engineer Skyrmion textures in Bose-Einstcin 
condensates. Analogous techniques could possibly also 
be used to couple FD atoms to em fields with the struc- 
ture of 3D monopole defects p^ . 

We have considered the em background field as a coher- 
ent classical field. This is consistent with the adiabatic 
approximation in the fermionic particle number fraction- 
ization, in which case the quantum fluctuations of the 
bosonic soliton are ignored H] . For a two-photon optical 



Raman coupling between the two species, the decoher- 
ence rate per atom is determined by the Rayleigh scat- 
tering rate ^ rmax(|7?,ip, |7?.2p)/A^, where F denotes 
the natural linewidth. This can be reduced by increas- 
ing the detuning A, provided that at the same time a 
sufficient laser intensity for the required tunneling rate 
Q, ^ TZiR.2/ Is. is available. In the case of a rf coupling 
the quantum effects of the em field can be safely ignored 
on the time scale of the experiments. 
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